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Abstract 

, We observed a slow relaxation of magnetoresistance in response to applied magnetic field in selectively doped p-GaAs-AlGaAs 
structures with partially filled upper Hubbard band. We have paid a special attention to exclude the effects related to temperature 
fluctuations. Though this effect is important, we have found that the general features of slow relaxation still persist. This behavior is 
interpreted as related to the properties of the Coulomb glass formed by charged centers with account of spin correlations, which are 
sensitive to an external magnetic field. Variation of the magnetic field changes numbers of impurity complexes of different types. 
As a result, it effects the shape and depth of the polaron gap formed at the states belonging to the percolation cluster responsible 
for the conductance. The suggested model explains both the qualitative behavior and the order of magnitude of the slowly relaxing 
magnetoresistance. 

PACS: 72.80.Ng 73.61 Jc 72.20.Ee 
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1. Introduction 

Recently JIJ] we reported observation of long time relaxation 
of (negative) magnetoresistance in p-type AlGaAs-GaAs quan- 
tum well structures where both wells and barriers were doped 
by Be. We argued that in these structures so-called A + centers 
- doubly occupied acceptors belonging to the upper Hubbard 
band - are formed in the well. The observed long time be- 
havior of magnetoresistance was explained by polaron effects 
involving the spin-correlated A + -centers J^. 

In the present paper, we report more detailed studies of slow 
relaxations in the same samples. The main point requiring more 
detailed measurements is that in our previous paper fl]] fluctua- 
tions of the sample temperature during measurements were out 
of proper control and account. Since temperature fluctuations 
may lead to a variation of the sample resistance our previous 
results and their interpretation need a proper verification. We 
are grateful to Z. Ovadyahy who attracted our attention to this 
problem. 

Investigation of role of the temperature fluctuations is the 
main goal of the present work. We will show that, although 
the temperature fluctuations indeed can produce a pronounced 
effect on the resistance, the observed long time relaxations of 
magnetoresistance have an independent source similar to that 
mechanism considered in (]]]. More careful analysis of the ex- 
perimental results has required, however, some modification of 
the theoretical model []]]. Namely, now we believe that two 
mechanisms can lead to a slow relaxation of the magnetoresis- 
tance. The fist one is the magnetic-field-induced shift of the 
chemical potential of the holes [Q]]. The second mechanism is 



due to a direct influence of the magnetic field on the aggregates 
responsible for formation of the polaron gap at the sites respon- 
sible for the conductance yfl. Contrary to yfl, we conclude that 
it is the second mechanism that explains the refined experimen- 
tal results. 

The paper is organized as follows. In Sec. [2] we briefly de- 
scribe the samples and experimental procedure and report the 
experimental results. These results are interpreted in Sec. [3] 
where the theoretical models are considered and compared with 
experiment. 

2. Samples, experimental procedure, and results 

We used the structures containing 10 GaAs quantum wells 
separated by Alo.3Gao.7As barriers. Thickness of both wells and 
barriers was 15 nm. Confining Alo.3Gao.7As layers had thick- 
ness of 20 nm. The growth procedure is described in detail 
in (]]]. Central regions with thickness of 5 nm of both wells and 
barriers were p-doped with Be (concentration 10 17 atoms/cm 3 ). 
The contacts were made by 2 min burning at 450° C in de- 
posited gold containing 3% of Zn. The samples were shaped 
as Hall bars. The resistance was determined from the voltage 
between the voltage probes at a fixed current of 1 nA. We stud- 
ied several samples cut from the same wafer. The samples were 
relatively low-Ohmic (10 s - 10 6 Ohms/D at 4 K) that is, in our 
opinion, due to the fact that the impurity band formed by A + 
centers is rather close to the valence band, the binding energy 
being ~ 10 meV [jj]. We have checked that for all our measure- 
ments the I - V curves were linear in the temperature domain 
4.2-1.35 K. The temperature dependence of the resistance is 
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compatible with the 2D Mott variable range hopping (see Fig. 1 
in EE), R k exp(r /r) I/3 , with T = 1000 K for the sample un- 
der investigation. The magnetoresistance curves were close to 
those reported in namely, the linear negative magnetoresis- 
tance crossed over to the quadratic positive magnetoresistance 
at stronger magnetic fields. 

Studies of the slow relaxations induced by controlled varia- 
tion of external magnetic field are specifically difficult because 
our thermometers are sensitive to magnetic field. Therefore, to 
investigate role of the temperature fluctuations we made special 
measurements of the sample resistance and its temperature in 
the absence of magnetic field. The corresponding curves are 
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Figure 1 : Evolution of the voltage across the sample and its temperature with 
time at H = The dashed curve corresponds to evolution of temperature while 
the solid curve - to evolution of voltage. The curves were drawn by the adjacent 
averaging of the experimental points (involving 10 points). 

presented in Fig. 1. One concludes that temperature fluctua- 
tions up to several mK indeed exist. We believe that these fluc- 
tuations are due to fluctuations of the pressure in the system 
pumped by a pre-evacuation pump without a sufficiently large 
damping volume, which should be at least by an order of mag- 
nitude larger than the Dewar volume. It is also seen that the 
resistance fluctuates keeping time with the temperature fluctua- 
tions - there is no lag between variations of the resistance and 
temperature. This fact allows us to conclude that at a given time 
the resistance of the sample is controlled by the temperature of 
the sample measured at the same instant of time. Accordingly, 
we can extract temperature derivative of the resistance and then 
estimate the corresponding fluctuation-induced contribution to 
the zero-field resistance at a known temperature as 



AR =a —AT , 
dT 



(1) 



From Fig. 1 it is seen that the variation of temperature by 1 mK 
leads to a variation of voltage and thus the variation of resis- 
tance is of the order of 0.2%. 

Shown in Fig. 2 is time dependence of magnetoresistance (or 
rather of the voltage between the probe contacts) obtained by 
ramping of magnetic field in time with very slow sweep rate 
(about 30 Oe/s ). Large circles show measured temperatures at 
H = 0, that correspond to 3 instants of time in Fig. [2] Note that 
the magnetoresistance curves of Fig. |2lwere close to those re- 
ported in Fig. 1 of our previous paper, yj] which were obtained 



at fast magnetic field sweep of ~ 200 Oe/s |4|]. Namely, the 
linear negative magnetoresistance crossed over to the quadratic 
positive magnetoresistance at stronger magnetic fields. As it 
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Figure 2: Magnetoresistance as a function of time. Small open squares repre- 
sent the voltage across the sample while the circles represent the temperatures 
at the instants of time when H = 0. Dashed lines connect the values of voltage 
and temperature at the times when H = 0. 

is seen from Fig. [2] the resistance at the successive instants of 
time when H — is different - it gradually decreases with time 
(by w 1%). At the same time, the temperature at that instants 
is almost the same showing initially weak decrease and then a 
weak increase within the interval not exceeding a 1 mK. These 
deviations of temperature would lead only to deviations of the 
resistance by « 0.2% (upward and downward). These mea- 
surements demonstrate that, in addition to thermal effects, the 
resistance tends to decrease after application of the magnetic 
field, the effect at the times ~ 1000 s being at the level of 1%. 

In Fig. 3, we present the behavior of the voltage across the 
probe contacts (and thus of the resistance) for two opposite 
current directions for the sample quenched during 20 min in 
a steady magnetic field of 0.7 T after subsequent switching off 
the field. In this case the temperature measurements are avail- 
able during all the time of the measurement. Thus we are able 
to correct the resistance curve for the temperature fluctuations 
at all experimental points. It is seen that the corrected voltage, 
U - {dU/dT)AT, slowly relaxes to higher values and is satu- 
rated at times ~ 800 s. It is important that the saturation value 
is equal to the resistance of the sample before application of the 
magnetic field. In other words, after application of the magnetic 
field the resistance slowly decreases with time but after switch- 
ing off of the field it is gradually restored to its initial value. It 
is also important that the restoration time is nearly equal to the 
time at which the sample was subjected to the external magnetic 
field. 



3. Discussion 

We believe that the results presented above evidence a spe- 
cific slow response of the sample resistance to applied magnetic 



2 




0,0- 



-0,1 -| — i — — i — | — i — | — i — | — i — | — i — | — i — i — i 

5800 6000 6200 6400 6600 6800 7000 7200 

t,S 

Figure 3: Evolution of sample resistance after switching off a steady magnetic 
field. The points correspond to values of the voltage for the two different current 
directions corrected for the effect of temperature while the curves are obtained 
by averaging procedure described above. The sharp peaks from the left corre- 
spond to the instant when the magnetic field is switched off. 

field. This response cannot be explained by temperature fluctu- 
ations, which would not lead to the observed contribution to the 
resistance, which monotonously depends on time. As it follows 
both from experiment and our estimates, a steady heating due 
to eddy currents induced in the Cu substrate by variations in the 
magnetic field is negligibly small. At the same time, the ob- 
served resistance gradually decreases after application of mag- 
netic field and tends to restore after switching it off. In addition, 
as it follows from direct measurements and our estimates, a pos- 
sible effect of temperature fluctuations would be about an order 
of magnitude less than the observed variation of resistance. 

3.1. Theoretical model 

Hole-impurity complexes. -. Let us fist review the previously 
suggested mechanism of slow relaxations HI, which is a gen- 
eralization of the one developed in yfl. We assume that some 
localized states (including A + -centers - acceptor atoms doubly 
occupied by holes and located in the wells) form bistable aggre- 
gates. The low energy states of these aggregates are almost de- 
generate, the transitions between the states can take place only 
due to many-electron processes that results in very long tran- 
sition times. Though the above aggregates do not belong to 
the percolation cluster, they can be polarized by the electrical 
charges located at the hopping sites. The polaron clouds formed 
by the aggregates for typical hopping sites lead to formation of 



a polaron gap at the Fermi level. As a result, conductance de- 
creases. 

An important feature of the material under consideration is 
that, in addition to doubly and singly occupied acceptors in the 
well (A + centers and A centers, respectively), there can also ex- 
ist neutral complexes consisting of negatively charged acceptor 
in the barrier and a valence-band hole located in the well, which 
is coupled to the acceptor due to Coulomb interaction. We will 
call such complexes the A centers. Analogof such complexes 
in a n-type conductor were introduced in J5[], and in what fol- 
lows we will exploit the scheme suggested in that paper. In the 
material under consideration, one has to consider acceptor com- 
plexes of 4 types: (i) positively charged A + -centers - doubly 
occupied (by holes) acceptors located in the wells and form- 
ing the upper Hubbard band; (ii) neutral A°-centers - singly 
occupied (by holes) acceptors located in the wells; (iii) nega- 
tive A~-centers - acceptors in the barriers; and (iv) neutral A°- 
complexes consisting of an acceptor in the barrier coupled to a 
valence-band hole in the well. 

Assuming that he dopant concentrations, N A , within the bar- 
rier and within the well are equal one has 

Na++N a o=N a -+N A o=N a . (2) 
- ' - ' 

well barrier 

Due to charge conservation we also have 

N A + = N A - . (3) 
In combination with Eq. (f2| we get 

N A o = N A o . (4) 

Since there are three unknown concentrations, N,, one needs 
one more relation in order to calculate the chemical potential, 
p. To derive this relation let us take into account the fact that 
a hole provided by an acceptor within the barrier is inevitably 
captured within the well. However, it can belong either to an 
A + or it can exist as a valence band hole coupled to a negative 
acceptor in the barrier forming an A center. In particular, a 
hole released by an acceptor in the barrier can be captured by 
an occupied acceptor in the well forming A + center provided 

e 2 

U A+ + = > U Aa . (5) 

K^ + {d n , + d b ) 2 /4 

Here <7, is the binding energy of the center of /th type, r is a 
distance between the two acceptors along the plane of the struc- 
ture, d w and d\, are the thicknesses of the well and the barrier, 
respectively. For simplicity we assume that kT <K U-, and that 
the dopants form (5-layers in the middle of the well and the bar- 
rier. The l.h.s. of Eq. ([5]) gives the energy gain due to formation 
of an A + -center (the second term describes the interaction of the 
A + -center with the A~-center within the barrier). One concludes 
that if Up < U A + all the acceptors within the well become the 
A + -centers (while the acceptors within the barriers all become 
the A " -centers). 
In contrast, if 

2e 2 

u # ~ Ua+ ~ 7T-T7\ (6) 
K(d w + d b ) 
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Figure 4: 



both A + -centers and A~-centers cannot be formed, and only A°- 
and A°-centers are present. Strictly speaking, in both of these 
limits the charge transport over the impurity band cannot take 
place. Such a transport requires partial "compensation" of the 
acceptors, which is possible provided 



< U/fi - U A * < 



2e 2 



K(d w + d h ) 



(7) 



In this case the relation between the concentrations of the A + - 
and A°-centers (and thus of the A°-centers) depends on the 
dopant concentration, Na, which controls the typical value of r 
in Eq. ©. Note that if Eq. © holds for r =* N A 1/2 the most part 
of acceptors within the well form A + centers while the num- 
ber of A centers is exponentially small (due to exponentially 



small probability to find the values of r > N A 
"compensation" takes place if 



-1/2 



< u 



k^N a 1 + (d w + d b ) 2 /4 



). The effective 



(8) 



In this case the relation ([5]) defines the average distance between 
A + centers and thus the relative densities of states of the upper 
Hubbard band (formed by A + - and A°-centers) and of the lower 
Hubbard band (formed by A - and A~ -centers), g u and gi, re- 
spectively. This distance, r corresponds to the equality relation 
in Eq. ©. 

Influence of magnetic field. -. The external magnetic field 
tending to align all the spins of the holes decreases the bind- 
ing energy of the A + -centers, Ua+, thus changing the balance 
between A + - and A°-centers. As it is known for A + -centers in 
GaAs/AlGaAs structures, the total spin of A + state is 2 while 
the spin of a separate hole is 3/2 (see JfJ). Thus the applica- 
tion of the magnetic field H leads to an effective increase of the 
Hubbard energy by the Zeeman energy term nsg\H\, 1±b being 
the Bohr magneton. As a result, the holes redistribute between 
the acceptors and the chemical potential is shifted as Q] 



Sfi ~ fi B g\H\ 



]u ~ gi 

gu + gl 



(9) 



Estimate for magnetoresistance. - . To compare the experi- 
mental results with our theoretical model let us reconstruct the 
estimate (jl for magnetoresistance. Consider a hopping site 
with equilibrium energy s (referred to the Fermi level) coupled 
to some two-level system (TLS) with inter-level splitting E in 
the equilibrium. The coupling potential U{R) (where R is a dis- 
tance between the hopping site and the TLS) is defined as the 
difference between the coupling energy for the two TLS states 
corresponding to the upper and lower levels of the TLS. If one 
creates an excitation at this hopping site (either an electron or a 
hole depending on the sign of e) the TLS splitting is changed as 
E -> E + U(R). If U is negative and \U\> E, the TLS changes 
its state with respect to the equilibrium one. Correspondingly, 
the excitation energy is also changed as s + U{R), and for an 
electron excitation its energy is lowered. For the hole excita- 
tion the same effect will take place if U > 0. As a result, the 
presence of a TLS leads to a formation of a polaron gap around 
the Fermi level having width of 2U and depth depending on the 
concentration of TLSs. The shape of the gap can be found as 
follows J^t]. Specifying the density of states for the TLSs with 
inter-level spacing E and relaxation time r as P(E,R, t) = P/t 
where P = const we find the distribution function of polaron 
shifts, U, from the equation 



T(U)dU 



InRdRP 
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U{R)P2nRdR In 



(10) 



This is just the probability to find a TLS providing the polaron 
shift between U and U + dU and located within the layer be- 
tween R and R + dR, which can switch between its states during 
the time of experiment, T max . Collecting contributions of all 
relevant TLSs, i.e., excluding all TLSs with polaron shifts < U, 
we obtain the shape of the polaron gap as 

F(U) = -27rP]n— f dU' U'R(U') [ ^- | . (11) 

Tmin Ju \dt/'/ 

Following our previous paper ^ let us consider many-electron 
aggregates(electronic TLSs) of chessboard type formed from 
the pairs with one occupied and one empty site. These building 
blocks can be regarded as pseudo-spins. The two states of the 
aggregate correspond then corresponds to opposite direction of 
all of the "spins". The transition between the two states can be 
due to either a coherent multi-electron hop, or due to motion of 
a "domain wall" separating parts of the aggregate with differ- 
ent phases of spin orientation. For this model, the size of the 
optimal aggregate turns out to be relatively large and its cou- 
pling to hopping sites is most effective when the hopping site 
is close to one of the sites belonging to the aggregate. In this 
case the coupling is just due to Coulomb interaction between 
the charges of the hopping site and the nearest site of the ag- 
gregate. The coupling to the rest sites of the aggregate turns 
out to be much weaker. In this case the distribution function P 
should be multiplied by the number of the sites belonging to 
the aggregate since the hopping site can be coupled to any site 
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within the aggregate. In what follows we will absorb this factor 
into P. For the Coulomb charge-charge interaction, U = e 2 /kR, 
the integrand in Eq. (fTTb is oc f/'~ 2 and, correspondingly, the in- 
tegral is controlled by its lower limit. The apparent divergence 
at U — » has a cut-off due to the fact that the small values of 
interaction energy U correspond to large distances between the 
hopping site and the aggregate. At large distance the Coulomb 
charge-charge coupling is replaced by a much weaker dipole 
coupling originating from a (random) dipole moment of the ag- 
gregate. Having this fact in mind we assume that the cut-off as 
the energy Sf, corresponding to the typical coupling energy be- 
tween the sites separated by the typical hopping length r/ t . The 
quantity £/, is just the width of the hopping energy band. Fol- 
lowing this reasoning we assume that the gaps with U < £/, do 
not effectively influence the hopping transport. As a result, the 
density of states near the Fermi level as a function of energy 
can be described by an interpolation relation as 



2nP 



2\ 2 

K / Ti 



•max (12) 



In the framework of the model the time scales T^^/min are 
controlled by the gate voltage protocol while In the present case 
case, the external perturbation is due to variation of the mag- 
netic field. An instant application of magnetic field leads to a 
shift of the chemical potential within the polaron gap. 

Let us estimate the corresponding change in conductance as- 
suming that 

S/i < fiagH < T « e h . 

From Eqs. ( TTZl i and (0 one can expect that energy shift of 
5e ~ 6p ~ fiBg\h\ will lead to a relative shift in the den- 
sity of states as well as in the i conductance of the order of 
(6fi/sh) 2 = (fiBgH/si,)) 2 . Therefore one can expect that the 
magnetoresistance in an oscillating magnetic field would form 
a pattern following the magnetic field. In addition to such a 
pattern, the experiments demonstrate a gradual increase of the 
conductance with time. 

To explain this increase we introduce an additional mecha- 
nism, which is due to a slow reconstruction of the polaron gap 
adjusting its depth to the magnetic-field-dependent number of 
the "active" sites forming the polaron cloud. The reconstruction 
can result from a direct influence of the magnetic field on the 
aggregates forming the polaron clouds. Indeed, the aggregates 
are formed from pairs of sites allowing electron (hole) transi- 
tions within the pairs. These transitions change types of the im- 
purity complexes as (A + , A ) -> (A , A + ), (A + , A") -> (A , A ), 
and (A , A ) — > (A~, A + ). Thus one concludes that the probabil- 
ity for the first configuration to be included into the aggregate 
is given by the product Na+Na" while for the second two con- 
figurations - by the product N A+ N 2 . In the region specified by 
Eq. (J8j, the magnetic field changes the relation between num- 
bers of the A + - and A°-centers. The effect of magnetic field on 
g u (related to Na+) and gi (related to A^o = Na<>) can be esti- 
mated from Eq. (0 for r - r c . 

Since the magnetic field changes the binding energy of an A + 
center, 6Ua* = P-BgH, it changes the characteristic distance, r c . 



Since Na* /Na — r 2 NA this change leads to a variation in the 
number of A + : 6Na+A + r 2 NA- Having this fact in mind and 
using Eq. (|5j one obtains: 



SN A * a 25U A +K[r 2 + (d w + d h ) 2 IA] 
N A * r 2 r e 2 



3/2 



5U A * 



Uso - Ua* 



As it was noted above, the concentration of the aggregates is 
controlled by the product Na*N a » = Na+(Na - Na*)- The varia- 
tion of this product, in its turn, is 



5N A+ (N A -2N A+ ). 



(13) 



Therefore one concludes that the effect of magnetic field de- 
pends on the sign of Na - 2Na* ■ If this sign is positive, that is 
the concentration of A + centers is relatively small, the concen- 
tration of the aggregates decreases with an application of mag- 
netic field (decreasing Na*)- In its turn, it leads to a suppression 
of the polaron effect and to a decrease of the resistance. If the 
concentration of A + -centers is large, the magnetic field can lead 
to an increase of the resistance. Since the experiment demon- 
strates a decrease of resistance as a result of the application of 
the magnetic field, we conclude that in our case the density of 
A°-centers exceeds the density of A + -centers. This assumption 
is supported by the fact that in our experiments the typical dis- 

— 1/2 

tances between acceptors, N A , were relatively large, at least 
N~ A V2 > (dy, +d b ) which is in favor of creating A centers rather 
than A + centers. 

The relative decrease of the density of states of the aggre- 
gates, 6P/P - SN A * /Na+ can be estimated from Eq. ( TT31 as 



SP ^ 2pg{\H\) 
P ~ U A o - U A * 



(14) 



where (\H\) is the time average of the absolute value of the mag- 
netic field. This ratio is equal to the relative (with respect to the 
total number of sites experiencing the polaron effect) number of 
the hopping sites where the polaron cloud is destroyed by the 
magnetic field. It is important that this effect is linear in (\H\) 
and can dominate over the quadratic effect mentioned above. 

For further estimates we have to specify the TLS distribu- 
tion function. A crude estimate for P in the case of electronic 
aggregates is HI 



Ne- 



(Np 2 )N l l 2 (e 2 lK P ) 



(15) 



Here is the number of pairs of sites forming the aggregate, 
the exponential is a statistical factor describing the probabil- 
ity to construct the bistable aggregate, p is the typical distance 
between the sites forming the aggregate. Here the first factor 
in the denominator describes the typical volume of the (2D) 
aggregate while the second - typical scatter of the TLS en- 
ergy splitting. The factor A depends on the competition be- 
tween the Coulomb interactions within the system and scatter 
of single-particle energies. Indeed, for weak Coulomb inter- 
actions the system is in its ground state and occupation of all 
single-particle states is given. It is the Coulomb correlations 
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that allow to have metastable configurations with close total en- 
ergies. Unfortunately, the number A is not known; it can de- 
pend on realization of the Coulomb glass. However one ex- 
pects that large relaxation times are available at not too large N 
and thus the exponential is not too small. Assuming N ~ 5, 
e~ AN as 0.1, p = « 100 nm (that is of the order of the 
typical hopping length) one estimates P as ~ 10 23 cm~ 2 erg -1 . 
Substituting this estimate in Eqs. dl2t . ( [14b and assuming that, 
Hg{H)/(U A „ - U A +) ~ 3 ■ 10~ 2 one obtains 



SG 
~G 



0.003 In 



(16) 



Here f max is the observation time while tn is given by the inverse 
sweep rate. This estimate by order of magnitude agrees with the 
experimental data. Note that the present estimate differs from 
that of [ 1 ] where the effect of magnetic field was attributed to 
partial suppression of the polaron gap by the magnetic field to 
an increase of the energy of the upper state of an aggregate. 
The suppression is due to influence of the magnetic field on the 
A + -centers making some configuration inaccessible even with 
account of the correlation energy. However, we underestimated 
a possible effect of an increase of the energy of the lower state 
involving A + centers, which can compensate the effect related 
to the increase of the energy of the upper level. In this way, we 
overestimated the slow-relaxing part of the conductance. On 
the other hand, the experimental result for this quantity was also 
overestimated because it was not corrected for the temperature 
fluctuations. We believe that the corrected experimental results 
and the suggested theoretical model are consistent. Namely, 
the model explains experimentally observed gradual increase 
of conductance with time. 

Thus we have discussed two different mechanisms of slow re- 
laxations (of magnetoresistance) in response to time-dependent 
external magnetic field, both related to the polaron effect. The 
first one is to some extent similar to the effect of the gate volt- 
age HIM 1 Illl2lll3ll . It is induced by magnetic -field-driven 
shift of the chemical potential. As it was noted above, this 
mechanism does not explain the observed gradual increase of 
conductance in time. We expect that this mechanism will be 
more pronounced and even dominant in the situation when the 
induced shift is larger than the width of the hopping energy 
band. 

The second effect is related to a direct dependence of the 
"active" sites responsible for the polaron effect on the magnetic 
field. According to our estimates, it is this effect that is respon- 
sible for the experimentally observed slow relaxing response. 

It worth noting that the slow relaxations in response to varia- 
tions of the gate voltage are usually not observed in doped crys- 
talline semiconductors - the experiments ]8|, ui llOj 111], [l2l 1 1311 
were performed using samples with significant amount of dis- 
order. We believe that there two reasons for that: (i) the TLS- 
induced polaron effects are much weaker that a direct influ- 
ence of the gate voltage on the DOS; (ii) relatively large sweep 
rates of the gate voltage probably lead to pronounced non- 
equilibrium behaviors. In the present experiments, the direct in- 
fluence of the magnetic field sweep is much weaker, and there- 
fore the effects induced by TLS polarons can be observed. 



To conclude, we observed a slow relaxation of magnetoresis- 
tance in response to applied magnetic field in selectively doped 
p-GaAs-AlGaAs structures with partially filled upper Hubbard 
band. We explain this behavior as related to the properties of 
the Coulomb glass formed by charged centers with account of 
spin correlations, which are sensitive to external magnetic field. 
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